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ABSTRACT 

2 ' Using stacks of Lya images of 2128 Lya emitters (LAEs) and 24 protocluster UV- 

^ \ selected galaxies (LBGs) at z = 3.1, we examine the surface brightness profiles of Lya 

^ ■ haloes around high-z galaxies as a function of environment and UV luminosity. We find 

that the slopes of the Lya radial profiles become fiatter as the Mpc-scale LAE surface 
densities increase, but they are almost independent of the central UV luminosities. The 
characteristic exponential scale lengths of the Lya haloes appear to be proportional 
to the square of the LAE surface densities (rLya cx ^lae)- Including the diffuse, 
^ ^ ^ extended Lya haloes, the rest-frame Lya equivalent width of the LAEs in the densest 

■ regions approaches EWq ~ 200 A, the maximum value expected for young (< 10^ yr) 
\ galaxies. This suggests that Lya photons formed via shock compression by gas outfiows 

or cooling radiation by gravitational gas infiows may partly contribute to illuminate 
the Lya haloes; however, most of their Lya luminosity can be explained by photo- 
ionisation by ionising photons or scattering of Lya photons produced in H II regions in 
and around the central galaxies. Regardless of the source of Lya photons, if the Lya 
haloes trace the overall gaseous structure following the dark matter distributions, it is 
not surprising that the Lya spatial extents depend more strongly on the surrounding 

■ Mpc-scale environment than on the activities of the central galaxies. 

^ \ Key words: galaxies: formation - cosmology: observations - cosmology: early uni- 

■ ■ ■ verse 
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1 INTRODUCTION 

It is believed that gaseous baryons around galaxies 
(circum-galactic medium or CGM) play an impor- 
tant role in determining the properties of galaxies 
via gas in fl ows / outflows across cosmic time (e.g., 
Keres et al.1 l2005l: iMori Umemural l2006l: IPekel et al.1 



20091: lOpDenheimer et aD l2010l : Ivan de Voort et al.1 l201ll : 
Fumagalli et al.l 12011 V Observationally, the CGM have 
been studied mainly using absorpt ion lines of back- 
ground QSOs or gal axies (e.g., Steidel et al.l I2OIOI : 

IProchaska et al.1 I2OIII : iRudie et all I2OI2I ). An alterna- 



^ Based on data collected at Subaru Telescope, which is operated 
by the National Astronomical Observatory of Japan, 
t E-mail: matsuda@astro.caltech.edu 



tive window of the CGM is Lya emitting halo, which 
should trace the rel atively dense gaseous structure 
aroun d a galaxy (e .g., Furlanetto et al. 2005: Goer dt et aP 



2OIOI: 

Rosdahl BlaizotI 



201^ . 



Faucher-Gigu erejt al. 



2010; 
biikstr 



Zheng et al] 



20111: 



Hultman Kramer! 



Difl[use, extended Lya haloes are likely a generic prop- 
erty of high-z galaxies. Based on a deep long- slit spec- 
troscopy of a QSO field, iRauch et all (|2008. ROSi ) detected 
diffuse Lya haloes by averaging the spectra of 27 Lya emit- 
ters (LAEs) z — 2.6 — 3.8, extending to a radi us of r ^ 
30kpc abo ve ^ 10~^^ ergs s~^ cm~^ arcsec"^. Steid el et al.l 
(|2011. siih detected more extended (r ~ 80 kpc) Lya haloes 
above a similar surface brightness limit by composites of the 
Lya images of 92 UV-selected galaxies in three protoclusters 
at 2; = 2.3, 2.8 and 3.1. Sll found that the Lya haloes can 
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be well fitted by an exponential profile and that the slopes of 
the Lya haloes depend weakly on the central Lya equivalent 
widths. 

What determines the Lya halo structure? The galax- 
ies in Sll are known to res ide in galaxy ove r- dense 
regions with ^gai ^ ^ " dsteidel et alJ |200Q| . l2005l : 
lErb. Bogosavlievic. Steidell l201lh . while most of R08's 
galaxies would reside in blank field. Moreover, the UV con- 
tinuum luminosities of the Sll's galaxy sample are > 10 
times brighter than those of the R08' s galaxies. To fill these 
gaps, we require observations covering larger dynamic ranges 
of environment and UV luminosity. Here, we examine the 
properties of Lya haloes around LAEs at z — 3.1 based on 
our wide-field deep Lya imaging data. 

In this paper, we use AB magnitudes and adopt cosmo- 
logical parameters, Qm = 0.3, Q^a — 0.7 and h — 0.7. In this 
cosmology, the Universe at z = 3.1 is 2.0 Gyr old and 1.0^^ 
corresponds to a physical length of 7.6 kpc at z = 3.1. 



2 OBSERVATIONS AND DATA REDUCTION 

We use the same data set as in iMatsuda et al ] (|201lh 
and lYamada et aD (|2012a. Y12l ). Full detail of the ob- 
servations and data reduction were described in Y12. 
Using a narrow-band fi lter, NB4Q 1 (CW=4977A and 
FWH M=77A. Havashino et ahlliool ) on Subaru Suprime- 
Cam (|Mivazaki et al.l I2OO2I I, we observed 12 pointings: 
Great Observatories Origins Deep Survey-North (GOODS- 
N), Subaru Deep Field (SDF), three fields in Subaru-XMM 
Deep Survey (SXDS-C, N, and S) and seven fields around 
SSA22 (SSA22-S bl-7). The SSA22 contai ns the protoclus- 
ter at z = 3.09 (jSteidel et al.lll998l . I2OOOI ). For the SSA22 
fields, we obtained broad-band {B and V) images in our 
observing runs. For GOOD S-N, we us ed archival raw B 
and IZ-band images (Capak et al.ll2004l ). For the SDF and 
SXDS fields, we used public, reduced ^ and V -band images 

(iKashikawa et al ] |2004l : lFurusawa et aLllioOsI ). 

We reduced the ra w data with sdfred fYa gi et aD 
I2OO2I : lOuchi et all HqqI) and IRAF. For photometric cal- 
ibration, we used t he photome t ric and s pectr ophotomet- 
ric standard stars (lOkd I199QI : iLandoltl (l992l ). We cor- 
rected the magnitudes using t h e Gal actic extinction map of 
ISchlegel Fi nkbeiner. David ("1998). We aligned the com- 
bined images and smooth with Gaussian kernels to match 
their seeing to a FWHM of 1.0^^ or 1.1^^ depending on the 
original seeing. We made BV images [BV = (2B + y)/3] for 
the continuum at the same effective wavelength as NB497 
and made NBc (continuum subtracted NB497) images for 
emission-line images. The total survey area after mask- 
ing low S/N regions and bright stars is 2.4 square de- 
grees and the survey volume is 1.8 x 10^ comovingMpc^. 
The 1-cr surface brightness limits of the NBc images are 
0.7- 1.2 X 10"^^ergs"^cm"2arcsec"^. 



3 RESULTS 

We use a sample of LAEs selected by Y12. The selection 
criteria are as follows. 

{1)NB497 < 25.73 (S/N > 6.6), and, 

{2a)BV - NB497 > 1.0 k B - Vc> 0.2 (K < 26.5), 
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Figure 1. Distribution of BV magnitudes and th e surface over- 
densit ies of 2 128 LAEs at z = 3.1 (red dots) from lYamada et "aTl 
(2012a . Y12D and 24 LBGs at 2 = 3.1 in the SSA22 protoclus- 
ter (blue stars) from Fsteidel et al.l |2003). The BV is the contin- 
uum magnitude at the rest-frame wavelength of ~ 1220A. The 
local over-density of LAEs (5lae = — ^)/^) is calculated by 
smoothing the LAE sky distributions with a Gaussian kernel of 
cr=1.5^ (or FWHM~ 5h~^ comoving Mpc). The dotted lines rep- 
resent boundaries between sub-samples for composites. 



or, 

{2h)BV - NBm > 1.3 {Vc > 26.5), 
where Vc is the line-corrected F-band magnitude. The B — 
Vc criteria is useful to further prevent the contamination of 
[O 11] emitters at z = 0.33 although we cannot apply it to 
the continuum faint sources, Vc > 26.5. Eor the robustness 
of the sample, we also apply an additional constraint that 
the observed BV — NB497 colour must be larger than the 4- 
a value. The previous spectroscopic observations show that 
the contamination of the foreground objects, mostly [Oil] 
emitters ai z = 0.33, is at most 1 ~ 3%, o r negligible for 
the LA E sample (jMatsuda et al.ll2005l . I2OO6I : lYamada etHI 
l2012bl ). 

From the total 2161 sample, we exclude 33 LAEs (1.5% 
of the whole sample) whose positions are within 20^^ from the 
edge of the images. In Figure[T] we show BV magnitudes and 
surface densities of the remaining 2128 LAEs. The BV mag- 
nitude represents the rest-frame UV continuum luminosity 
near 1220 A. The surface over-density (^lae = (S — S)/S) 
is derived by smoothing the LAE sky distribution with 
a Gaussian kernel with a size of a = 1.5^ (or FWHM 
^ bh~^ comoving Mpc) as used in Y12, whose scale is close 
to the median distance between the nearest neighbours of 
LAEs in the blank fields (SXDS, GOODS-N and SDF). We 
confirmed that the result does not change significantly if we 
use slightly different kernel sizes. We divide the LAE sample 
into four sub-samples based on the LAE surface density or 
the central BV magnitude (see Figure [T] and Table [T]). We 
chose the boundaries between the sub-samples as follows. 
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Table 1. Properties of composite sub-samples 



Sample 


N 


"lae 


NB^, 
(mag) 


BV^ 
(mag) 


NB- 
(mag) 


BV 
(mag) 


(A) 


(A) 


(10-18) 


NBc 

(kpc) 


^BV 
(10-18) 


'-BV 

(kpc) 


2.5 < S^^^ < 5.5 
1.5 < S^^^ < 2.5 
0.5 < < 1.5 
-1 < < 0-5 


130 
273 
861 
864 


3.5 
1.9 
0.9 
0.1 


25.50 
25.36 
25.40 
25.35 


27.17 
26.94 
26.95 
26.97 


24.26 
24.49 
24.70 
24.77 


26.73 
26.47 
26.64 
26.68 


88 
81 
79 
84 


183 
117 
112 
109 


0.7 
1.4 
1.4 
1.5 


20.4 
13.2 
10.7 
9.1 


0.3 


10.3 


21 < BV^ < 25 

25 < BV^ < 26 

26 < BV^ < 27 

27 < BV^ < 29 


31 
203 
894 
1000 


0.5 
0.8 
0.9 
0.9 


23.54 
24.58 
25.33 
25.53 


24.59 
25.66 
26.60 
27.59 


22.39 
23.78 
24.65 
25.00 


23.91 
25.05 
26.15 
27.39 


50 
51 
61 
126 


76 
61 
75 
170 


3.2 
1.8 
1.6 
0.7 


14.3 
13.2 
10.2 
11.9 


0.9 
0.8 
0.3 


18.8 
13.3 
14.0 


Protocluster LBG^ 


24 


3.7 


26.19 


25.31 


23.05 


24.18 


8 


53 


3.2 


18.2 


5.0 


7.5 



"The LAE (surface) over-densities calculated by smoothing the LAE spatial distributions with a Gaussian kernel of cr=1.5^ (or 
FWHM~ 5/i-i comovingMpc). 

^The magnitudes (AB) measured with 2" diameter aperture photometry. 

*^The magnitudes (AB) measured with isophotal apertures determined on the NBc (Lya) images with a threshold of 
31.5 ABmagarcsec-^ (or 8 X 10-^° erg s-i cm-^ arcsec-^). 

'^The rest-frame Lyo; equivalent widths from the 2" diameter aperture photometry. 
^The rest-frame Lya equivalent widths from the isophotal aperture photometry. 

•^Best-fit parameters assuming a surface brightness profile of S{r) = Cnexp{—r/rn), where Cn is a constant in units of 10-18 erg 
s-i cm-^ arcsec-^ and is a scale length. 

^LBG sample in the SSA22 protocluster from lSteidel et al.l (|2003h . 




Figure 2. Composite images of z = 3.1 LAE sub-samples and protocluster LBG sample. The upper panels are Lyo; images and the 
lower panels are BV continuum images. All the stacked Lyo; images show extended Lyo; haloes. The size of the images is 20^^ x 20^^ 
(~ 150 X 150 kpc^). The white horizontal bars indicate the angular scale of 100 kpc (physical) at z = 3.1. The images are displayed with 
a logarithmic scale with a range between 26 and 33 ABmag arcsec-^. 



We first divide the sample into two sub-samples near the 
peaks of the density or BV magnitude distributions. Then 
the dynamic ranges of the high density or BV bright samples 
are much larger than those of the low density or BV faint 
samples. Therefore we further divide the high density or BV 
bright samples so that sub-samples have basically a dynamic 
range of ABV=1 or A(5lae=1- However, as the densest or 
BV brightest samples have too small numbers of sources, we 
use larger dynamic ranges {ABV^A or A5lae=3) so that 
they have enough number of sources. For comparison, we 
also use a sample of 24 LBGs at z = 3.06 — 3.13 (i.e. Lya 
reds hift ran^e cove r ed by NB497) in the SSA22 protocluster 
from lSteidel et al.l (|2003l ). 

Figure [2] shows composite NBc (Ly<^) and BV (Contin- 
uum) images of the LAE sub-samples and the protocluster 
LBG sample. We make the composite images as follows. We 
cut out 40'' X 40'' (or 300 x 300 kpc at z = 3.1) images cen- 
tred at the sources. We then stack these images with me- 
dian to erase unrelated foreground or background sources 
near the sources without masking before stacking. Median 



should also be more robust than average if th e samples in- 
clude rare luminous Lya blobs (LABs, e.g., ISteidel et alj 
|200Q| V The 1-a surface brightness limits, estimated from the 
fluctuations of the outer part of the composite images (> 10" 
away from the source), reach ^ 30.5 — 32.5 magarcsec"^ 
or ^ 0.3 — 2 X 10~^^ ergs s~^ cm~^ arcsec"^, depending 
on the number of stacked sub-samples. These depths are 
equivalent to those achieved with ^ 30 — 1000 nights ob- 
servation with 8-m class telescopes and is > 10^ times 
fainter th an the sky brightn ess at this wavelength on Mauna 
kea (e.g., iKrisciunas et al.lll987l V All the stacked Lya im- 
ages show significant extended Lya haloes. The sub-samples 
with brighter UV continuum show extended UV continuum 
haloe£]. For the LAE sub-samples, aperture corrections for 
NBc from the 2" diameter aperture to isophotal aperture 
photometry (above 31.5 ABmag arcsec"^ on NBc images) 

1 Similar diffuse (and/or clumpy) continuum s t ructure is also 
seen in giant LABs (e.g., iMatsuda et al.l [20071 : iPrescott et al.l 

[mil). 
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Figure 3. Lya and continuum surface brightness radial profiles of stacked LAEs and LBGs. Plots (a) and (b) are Lyo; profiles of 
surface density and BV magnitudes divided sub-samples respectively. Plots (c) and (d) are UV continuum profiles of surface density and 
BV magnitudes divided sub-samples respectively. The grey shades show the range of point spread function (PSF) of the 12 pointings, 
normalized at the smallest bin of lowest density or faintest UV luminosity sub-samples. All the Lya haloes are much larger than the PSF. 
For continuum, it is difficult to examine the halo profiles due to their faintness except for the sub-samples with brighter UV continuum. 
The dashed lines are fitted exponential profiles, S{r) = Cnexp{—r /rn)^ where Cn is a constant and rn is a scale length. In order to avoid 
effect from the PSF and central bright cores, only bins larger than r = 3'' are used for the profile fitting. The error bar shows a 1-a 
uncertainty for each bin. 



are ^ 0.5 — 1.2 mag, while those for BV are ^ 0.2 — 0.7 mag. 
We have confirmed that the results do not change signifi- 
cantly if we use different stacking methods (e.g., 2 or 3-a 
clipping average). 

The composite Lya and continuum surface brightness 
profiles are shown in Figure [3l We also plot the range of 
the point spread function (PSF) of the 12 pointings, nor- 
malized with the data at the smallest radius for lowest LAE 



surface density or faintest UV luminosity sub-samples. The 
1-a uncertainty for each bin is estimated from dividing the 
1-a surface brightness limit of each composite image with 
^/]4, where A is an area of a bin annulus in unit of arcsec^ 
(i.e., assuming Poisson statistics). To quantify the shape 
of the haloes, we fitted them using an exponential profile, 
S{r) — Cnexp{—r/rn), where Cn is a constant and rn is a 
scale length, by following Sll. In order to avoid contribu- 
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Figure 4. Lyo; scale length as a function of (a) LAE surface density and (b) central UV luminosity based on BV magnitude measured 
with 2" diameter aperture photometry. The symbols are same as in FigureO The dashed lines are r = 0.6 x (1 + ^lae)^ + 8.4 kpc (left) 
and r = 12.4 kpc (right). These trends suggest that the spatial extents of the Lya haloes are determined by the surrounding Mpc-scale 
environment rather than the central UV luminosities. The upside down triangles are not main sub-samples shown in Tabled but parts 
of sub-samples for checking these trends (see text). The error bars represent 1-cr uncertainties of the exponential pr ofile fitting shown in 
Figure O The star-formation rate (SFR) is converted from the BV magnitude without dust attenuation correction (|Kennicuttlll998l V 



tion from the central galaxies and effect from the PSF, we 
use only data larger than a radius of 3'' for the fitting. The 
best-fit parameters are shown in Table [1] 

We check consistency between our work and the Sll's 
work using the protocluster LBG sample. The Lya scale 
length of the protocluster LBG sample (r = 18.2 ± 2.2 kpc) 
is similar to that of the Sll's full sample (r = 17.5 kpc), both 
of which are derived from median stacking. In contrast, the 
UV-continuum scale length of the protocluster LBG sam- 
ple (r = 7.5 ± 2.0 kpc) seems to be larger than Sll's one 
(r = 3.4 kpc). This difference would come from the different 
radius ranges used for the profile fittings. We use data larger 
than a radius of 3^^ for the fitting of the haloes. However, the 
surface brightness of the UV-continuum profile with r > 3^^ 
are fainter than ^ 31 ABmag arcsec^ and almost out of the 
range studied in Sll (as shown in their Figure 7). If we 
use data within r < 3'' or above 31 ABmag arcsec^ , the UV 
continuum scale length of the protocluster LBG sample de- 
creases to r 4 kpc and approaches the Sll's result. 

In Figure m we plot the characteristic exponential scale 
lengths as a function of the LAE surface over-density and 
central BV magnitude. The Lya scale lengths appear to 
be roughly proportional to the square of the LAE surface 
density, r oc (1 + ^lae)^ oc Slae while there is no clear 
trend with the central BV magnitude. We check that the 
LAE sample in a peak over-density with ^lae = 4.5 — 5.5 



has a Lya scale length of rLya = 29.2 ± 4.8 kpc and also 
follows this trend (orange upside down triangle in the left 
panel). The Lya scale length of the BV brightest {BV = 
21 — 25ABmag) sub-sample seems to have slightly higher 
value than other BV magnitude divided sub-samples. After 
excluding 6 sources in BV — 21 — 24ABm a^, which may 
contain QSOs or AGNs (jSteidel et al.ll2002l \ the LAE sub- 
sample with BV = 24 — 25 ABmag has a Lya scale length 
of TLycK = 12.0 ±0.7 kpc (orange upside down triangle in the 
right panel), and is consistent with the average Lya scale 
length of the 5U-magnitude divided sub-samples. The Lya 
scale length of the protocluster LBGs is close to that of the 
LAE sub-sample with a similar over-density (left panel) but 
larger than those of the BV magnitude divided sub-samples 
(right panel). These results suggest that Lya scale length 
is a strong function of the Mpc-scale LAE surface density 
rather than the central UV luminosity. 



4 SUMMARY AND DISCUSSION 

Using samples of 2128 LAEs and 24 protocluster LBGs at 
z = 3.1, we examine Lya surface brightness profiles as a 
function of environment or central UV luminosity. Our re- 
sults show that the Lya scale lengths are roughly propor- 
tional to the square of the Mpc-scale LAE surface densities 
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(^Lya oc Slae)i t)ut are almost independent of the central 
UV luminosities. 

The apparent profiles of the Lya haloes should be deter- 
mined by a combination of gas distribution around galaxies 
and source(s) of Lya photons (e.g., R08; Sll). Including the 
Lya haloes, the rest-frame Lya equivalent widths of all the 
density divided LAE sub-samples are larger than EWq ^ 
100 A expected for star-forming galaxies with continuous 
star formation lasting > 10^ yrs (e.g., Sll) and that in 
the over-dense region almost reaches the maximum value of 
^ 20 A for very young (< lO'^yrs) galaxies (Chariot k Faj] 
1 19931 ). Therefore, it is possible that Lyg photons formed via 
shock compression by gas outflows (jMori Umemu'rall 20061 : 
lYamada et al.l 2012bll or cooling radiation by gravitational 
gas inflows j Coerdt e t al.|[2Qlol : 'Paucher-Gi guere et al.ll2010l : 
iRosdahl k Blaizot 20n|)" partly contribute to illuminate the 
Lya haloes especially in over-dense environment. However, 
more than half of the Lya luminosity from the haloes can be 
explained by photo-ionisation by ionising photons or scat- 
tering of Lya photons produced from Hll regions in and 
around the central galaxies in any case. 

Although the CGM structur e could be modified by 
gas inflow / outflow processes (|Mori Umemural 120061 : 
iGoerdt et all I2OIOI I. the overall struct ure is likely to fol- 
low the dark matt er distributions fe.g. jYauma et aDl201ll : 
iRudie et al.1 I2OI2II . From cosmological numerical simula- 
tions, ^Zh^rig~eriiL| (| 20 111 ) predicted that Lya surface bright- 
ness profiles are made from a combination of dark matter's 
one-halo and two-halo contributions with a transition scale 
of r 50kpc for a ^ 10^^ M© halo. Accounting for this 
model, our results can be qualitatively interpreted as fol- 
lows: the observed steeper slope of the Lya haloes in under- 
dense regions mainly represent the one-halo term while the 
flatter slope in over-dense regions is dominated by the two- 
halo term. In this case, it is not surprising if the Lya spatial 
extents are dete rmined by the surrou nding Mpc-scale en- 
vironment fe.g., lAbbas Shethll2005l V This interpretation 
seems to be supported by the known tight links between 
individually detected LABs and the surrounding large-scale 
structure: giant LABs have been preferentially discovered in 
galaxy over-dense regions on Mpc-sca l es (e.g. ISteidel et al.l 
'2000; Matsuda et al. 2004. 2009*. '201 1': 'Prescott et al.l l20o'8[ 
lYang et al. 2010 : Erb. Bogosavlievic. Steidel 20111 

This study is the first attempt at investigating Lya 
haloes around high-^ galaxies in a wide range of environ- 
ment and central UV luminosity. However, we lack a sample 
of UV-bright galaxies with deep Lya imaging in average to 
moderate galaxy density environment as shown in Figure (TJ 
It would be interesting to test if the Lya haloes around 
UV-selected galaxies also show such strong environmental 
dependence. In addition to Lya haloes, our data hint ex- 
tended continuum haloes at least for UV continuum bright 
LAE sub-samples. Although we focus on only Lya haloes in 
this paper, further studies of the UV continuum haloes will 
be important to investigate Hii regions in the CGM. More- 
over, future ^ 30deg^ Lya imaging survey with Subaru / 
Hyper Suprime-Cam will provide us with 10 k LAE sam- 
ple (per redshift slice) at z = 2 — 7 and enable us to study 
the redshift evolution of the CGM structure via Lya haloes. 
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